arXiv:1507.07383vl [astro-ph.EP] 27 Jul 2015 


Equilibrium rotation of semiliquid exoplanets and satellites 

Valeri V. Makarov 

vvm@usno.navy.mil 

US Naval Observatory, 3450 Massachusetts Ave NW, Washington DC 20392-5420, USA 
Accepted . Received ; in original form 

ABSTRACT 

A wide range of exoplanet and exomoon models are characterized by a finite 
average rigidity and a viscosity much lower than the typical values for terrestrials. 
Such semiliquid bodies may or may not have rigid crusts with permanent figures. 
Unlike planets with solid mantles and Earth-like rheology, semiliquid bodies can 
be captured into stable pseudosynchronous spin resonance, where the average rate 
of rotation is higher than the synchronous 1:1 resonance. Two basic conditions 
are derived for capture of planets with a triaxial figure into pseudosynchronous 
rotation, one related to the characteristic tidal wave number (the product of the 
tidal frequency by the Maxwell time), and the other to the orbital eccentricity. 

If a semiliquid object does not satisfy either of the two conditions, it is captured 
into the synchronous resonance. For nearly axially symmetric bodies, only the first 
condition is in place, and the other is much relaxed, so they should predominantly 
be pseudosynchronous. It is also pointed out that the equilibrium pseudosychronous 
rotation rate can not reach the widely used asymptotic value from the constant time 
lag model but is in reality closer to the synchronous spin. 


Introduction 


Most of the Solar system bodies subject to significant tidal interaction are captured into 
the 1:1 spin-orbit, resonance, also known as synchronous rotation, when the body’s average 
sidereal spin rate is exactly equal to the orbital motion. As a result of this resonance, the 
body is permanently facing the perturber with the same side, a phenomenon well known to us 
from the Moon. A notable exception to this rule is Mercury, which is captured in a stable 3:2 
resonance, whi ch makes one Mercurian solar day as long as two Mercurian solar years. As was 
shown by IHut (119801 1. synchronous rotation and a circular orbit is the terminal point of the 
tidal spin-orbit evolution; however, higher stable spin-orbit resonances a re possib le and even 
likely for super-Earth exoplanets of presumably terrestrial composition (IMakarov fe Berghea 
2 0141 ). Like Mercury, such planets, once caught in a super-synchronous resonance, should stay 
there indefinitely long, unless an external impulse, such as a giant impact, drives the planet 
out of this locally stable state. Most of the Solar system moons, on the other hand, are locked 
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in the 1:1 resonance, and no credible evidence was found of an equilibrium rotation with a 
non-resonant rate above the synchronous rate. 


Investigation of equilibrium near-synchronous rotation states was pioneered by Goldreich 


(1196611 for the ad hoc constant-Q tidal model, in which the tidal torque is independent in 
magnitude of the perturbation frequency. For a long time the work-horse of the tidal dynamics, 
the Constant Time Lag (CTL) model of tidal dissipation assumes that the secular tidal torque 
acting on the perturbed body can be approximated by a linear dependence on tidal frequency 
v with a negative slope ( analogous to a regular f orce of friction), but shifted upward for a 
non-vanishing eccentricity (IMurrav fe Dermottlll999f) . As a result of the shift, the secular tidal 
torque is positive (accelerating in the prograde sense) at the point of resonance, i.e., where the 
spin rate exactly equals the mean orb ital motion, 9 = n. The secular term, which acts as friction 


in the balance of forces, vani shes at (jHutlll98ll : iRodrigues et al .1120081 : iFerraz-Mello et al.ll2008 
Williams fe Efroimskv 2012 1 

6 ~ (1 + 6 e 2 )n (1) 


If the CTL model were applicable, in the absence of other forces, the moons of the Solar 
system would have been involved in a slow prograde rotation relative to their host planets. 
The commonly accepted explanation why this does not happeijj] refers to the presence of a 
geometrical triaxial torque caused by the elongation of the permanent figure. This torque 
emerges whenever the axis of the smallest inertia moment begins to deviate from the instanta¬ 
neous direction to the perturber; thus, it acts as a restorative force directed toward the point 
of exact synchronism. 


Pseud o-sy nchronous equilibrium emerges in other, more advanced and complex tidal mod¬ 
els as well. Ferraz-M cllo ( 2013 1 considered a rheological model dominated by Newtonian creep 
mechanism of dissipation and found an equilibrium tide at a super-synchronous rate of ro¬ 
tation, although the shift turned out to be frequency-dependent in that case, as opposed to 
simpl y 6 e 2 n. The hybrid viscoelastic model with self-gravitation and Andrade creep developed 
by Efr oi msk v (12012a .bl). also permits super-synchronous equilibria in principle. However, this 
state was shown to be practically impossible for bodies of t errestr ial com positio n, character¬ 
ized by large effective vi scosities and long Maxwell times dMakarov fe Efroimskvl 120131 1. see 
also (ICorreia et al. 20L4). On the other hand, a wide range of plausible rheological models 
for exoplanets and exomoons, as well as for smaller satellites of the Solar system, can be best 
represented by the viscoelastic paradigm with low values of viscosity but finite rigidity. The 
first important class of such objects is a terrestrial body that has been heated up by strong 
tidal forces past the point of solidus in the deeper layers of the mantle. 


1 There are lingering doubts about the exact synchronism of Europa though ( Geissler et al. 1 
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2. Torques 


In the following, we consider tides and forces acting on one of the components in a system 
of two orbiting bodies. The component under consideration is called in this paper “perturbed 
body”, or “rotating body”, and in practice can be a planet or a satellite. The other component 
is called “perturbing body”, or simply“perturber”, and it can be a star or a planet. The 
possible tides or deformations of the perturber are not involved in the analysis; it is effectively 
a point-mass with only a central gravitational action on the perturbed body. 


Comprehensive equations for the polar tidal torque (i.e., the torque component directed 
alo ng the axis o f rotation), including the fast oscillating terms, can be found, for example, 
in (IMakarov fe Bergheal 120141) . Here we only summarize the main equations necessary for 
the following analysis of equilibrium rotation. The non-polar component of perturbation are 
ignored in this analysis, which is equivalent to setting the orbit obliquity on the perturbed 
body’s equator to zero. The instantaneous torque acting on a rotating body is the sum of 
the triaxial torque caused by the quadrupole inertial momentum, and the tidal torque caused 
by the dynamic deformation of its body. Planets and moons of terrestrial composition have a 
permanent figure, which results in the three principal moments of inertia A, B , and C (the C 
denoting traditionally the largest moment of inertia) being unequal to one another. The degree 
of elongation of the permanent figure, or the departure from the dynamical axial symmetry, 
can be represented by the parameter (B — A)/C, which ranges from « 0.1 for smaller moons to 
pa ICE 5 for the larger rocky planets. The angular acceleration is related to the polar components 
of these torques by 

a ^tri + Ttide /n ^ 

6 = WP (2) 

with 6 being the rotation angle of the perturbed body reckoned from the axis of its largest 
elongation, corresponding to the moment A, M 2 is its mass, and £ being the coefficient of 
inertia, C/(M 2 R 2 ). Possible variations of the moments of inertia on longer time scales are 
not considered in this paper. The equations of the triaxial torque Ttrt have been exploited 


in th e literature on many occasions flDanbvlll962l; iGoldreich fe Pealall968t iMurrav fe Dermott 


19991) . Neglecting the precession and nutation of the planet, the angular acceleration due to 


the triaxial torque is 


#tri — —~ 


3 B 
2 


A <r 


n 2 sin 2(9 — /), 


( 3 ) 


C r 3 

where r is the instantaneous distance of the perturbed body from the perturber, a is the semi¬ 
major axis, n is the orbital mean motion, and / is the instantaneous true anomaly. This formula 
can be more conveniently rewritten as an infinite series of mean anomaly (A4) harmonics: 


$tri — — 


3 B-A 


C 


-n 


G 20 O- 2 ) sin(2 6 — jA4), 
j 


( 4 ) 


where Gi pq are Kaula’s functions of eccentricity related to Hansen’s coefficients through G 200 - 2 ) = 
_ 0 2 

Xj ' . Note that the summation in Eq. [4] is over j = — 00 ,..., + 00 , and £t 2 o(- 2 ) = 0. 
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A simplifi ed e quation for the angular acceleration caused by tidal torque can be written 
as f Efroimskv 2012 a. Eqs. 106-108) 


6. 


3 Mi 


TIDE — 


2 £M 2 


-n 


— ) ^2 ^ 2o <?( e ) ^2 < ^ 2 ° j ( e ) 

9 =— i i=—i 


[A' c (2, i/ 220 q) Sign(cu 220g ) cos ((g - j)7W) - K a ( 2, i/ 220g ) sin ((g - j)A4)], (5) 

assuming M 2 <C Mi where M\ is the mass of the perturbing body, and 


A c (2, u 220q ) = fc 2 (^ 22 o g ) sin e 2 (^ 22 o g ) , 

K s (2, V220q) = k 2 (v 22 o q ) COSe 2 (z/ 22 0 q) , (6) 

fc 2 being being the degree-2 dynamical Love number, e 2 (z/ 22 o 9 ) the positive definite, mode- 
dependent tidal phase lag and u 220q being the positive definite physical forcing frequency related 
to the semidiurnal Fourier mode through 


^220 q — I w 220g| — | (2 + q)n — 29\. 


(7) 


One of the functions, K c , has been scrutinized in the literature flEfroimskvi l2012a .bl). Histor¬ 
ically, it is denoted as k 2 /Q, also dubbed kvalitet elsewhere, the Q value being related to the 
amount of energy dissipated over one cycle of rotation. Secular (i.e., non-averaging over one 
orbital period) acceleration emerges from the K c term when q — j. 

Within the simplifying assumption of a spherical, incompressible, and homogeneous body, 
and a pure ly analytical approach ad c tpted in this pap er, the corresponding kvalitet can be 
written as (lEfroimskv fe Lainevl 120071; lEfroimskvll2012alJbl ) 


K c ( 2,u) 


3 A 2 vT M ^s(y) 

2 (3R(z/) + A 2 vt m ) 2 + A 2 (V) 


where 


A 2 


Q being the gravitational constant, and 


, 19 R A n 

An -— 

6 QMl 


/ Q/'jr \ 

Sft(z/) = ur M + A- Q r M TA a COS J r(l + a) 
= -1“ ^- a r M rX a sin (^) T(1 + a) 


( 8 ) 

(9) 


( 10 ) 


The 3?(z/) and A(z/) dimensionless functions are the real and imaginary parts of the normalized 
dimensionless complex compliance, respectively, multiplied by tidal wave number vtm to avoid 
numerical singularity at zero frequency. In these equations, R is the radius of the perturbed 
body, M 2 is its mass, tm is the Maxwell time, ta is the Andrade time, /i is the unrelaxed shear 
modulus. The frequency v is used hereafter as a shorthand version of z/ 220( j, so the functions 
remain dependent on both the frequency and the tidal mode. The dimensionless Andrade 
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parameter a takes values from ~ 0.16 for a partial melt to ~ 0.4 for cold minerals. The 
parameter A, sometimes call ed effect ive rigidity in the literature, is related to the definition of 
the static Love number (jGoldreich fc Pealei ll96 8l ). The second terms of the complex compliance 
components represent the contribution of the Andrade creep mechanism. 


These basic equations greatly simplify for the case of semiliquid bodies considered in this 
paper. Recalling that the Maxwell time is defined as tm = 77 //x, where 77 is the effective 
viscosity, and that the Andrade creep is perhaps insignificant for such bodies, one can use the 
constraints ta = 00 and tmu <C 1. These simplify the secular kvalitet to 

K cM = k 2 {v) sin e 2 (u) = 2 ( 1h ^) 2 + 1 > ( n ) 


where A = TmV is a dimensionless quantity hereafter termed the tidal wave number. This 
function is strongly dependent on tidal frequency in narrow vicinities of spin-orbit resonances 
(2 + q)n = 26 for integer q. The characteristic kink-shape of the near-resonant torque has the 
largest amplitude in the 1:1 resonance for small and moderate orbital eccentricities. A steep 
decline between the two peaks (Fig. Q]) occupies a narrow band of frequencies for realistic 
rheologies. Despite the relatively small amplitude of the kink (compared to a typical amplitude 
of the triaxiality-caused torque), it acts as an efficient trap for a planet trying to traverse the 
resonance. The maximum kvalitet is reached at 


where its value is 


^max 

1 

(12) 

i + m 2 ’ 

-^"c(2, ^max) 

_ 3 A 2 
~ 41 +A 2 ' 

(13) 


3. First condition of pseudosynchronous equilibrium 


As is evident from Eq. [12l the locale of maximum tidal force in frequency is inversely 
proportional to the effective Maxwell time. Relatively cold terrestrial mantles are characterized 


by large values of viscosity (e.g., iDoremusI 120021 ) . and, thus, long Ma xwell times. For exa mple. 


the tidal response of Mercury was modeled with ar« = 500 years (iNovelles et al. 2014), and 
the suggested value for the Moon is ~ 10 years as estimated from the observed tilt of the 
mean principal axis flMakarov 2013 1. Icy satellites, on the other hand, can have much shorter 
Maxwell times owing to the relatively low viscosity of pressurized ice ( rs-/ 10 15 Pa- s). The 
slope of the central part of the secular torque between the opposing peaks (Fig. |T] left), which 
can be roughly approximated with an inclined line, is there fore proportional to the Maxwell 
time. Within this segment of the kink-shape, the secular torque acts as a quasi-linear friction 
force on any deviations of the spin rate from the equilibrium value. The point of equilibrium, 
or zero secular torque, is associated with the root of the torque curve, which is offset from 
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rate of rotation 91 n-1 


rate of rotation 91 n-1 


Fig. 1.— Angular acceleration of Titan caused by the secular component of tidal torque (Eq. 
0 with the assumed parameters listed in Table [Hand with a Maxwell time of 0.1 days. Left: a 
large-scale view of the dependence showing the main “kink” corresponding to the 1:1 resonance 
term, c^oo- Right: the same dependence on much smaller scale to demonstrate the offset of 
null torque from the exact 1:1 resonance. The kinks associated with the other resonances 
9/n — 0.5,1.5, 2.0, 2.5,... are too small in amplitude and smoothed out to be discerned. 

the point of synchronous resonance because of the positive, nearly constant, component due 
to the second-largest term q — 1 for a finite eccentricity (Fig. [T] right). The tidal torque 
drives the rotating body away from the point of synchronism 1^2200 — 0 toward the point of 
equilibrium z / 2 200 = ^pseudo when the frequency takes values within the narrow segment between 
these points. Obviously, the point of pseudosynchronous equilibrium, where the secular torque 
vanishes, can not be outside of the linear segment bracketed by the opposite peaks. Therefore, 
the asymptotic value of (1 + 6e 2 )n can only be realized when this frequency is much smaller 
than the location of the peak r^ 1 (1 + M 2 ) -1 , he., in liquid or semiliquid bodies. For inviscid 
planets with Maxwell time comparable to, or greater than, the orbital period of the perturber, 
the point of equilibrium collapses to very small values. This practically eliminates the feasibility 
of pseudosynchronous equilibrium, as any small perturbation, including physical libration, will 
drive the body into the globally stable synchronous rotation. 

Fig. [ 2 ] shows the relative, normalized frequency of pseudosynchronous equilibrium as a 
function of tidal wave number for the Titan model (Tabled]). This graph is obtained by com¬ 
puting the roots of Eq. [5] for the secular part of tidal acceleration (q — j) and for different 
values of tm- The previously known point of (1 + 6e 2 )n from the CTL formalism, indicated 
by the horizontal line in the graph, is only an asymptotic limit achieved at A —> 0. Realis¬ 
tic viscoelastic bodies with self-gravitation attain the pseudosynchronous equilibrium at any 
rotation rate between n and (1 + 6 e 2 ) n depending on the effective wave number tm n. 

As is shown in Fig. [2] the point of equilibrium becomes close to the CTL limit at wave 
numbers as small as 0.1. This implies a Maxwell time of 10 days or shorter for the Titan model 
(Tabled])' The fact that Titan is in synchronous rather than pseudosynchronous rotation can 








equilibrium tide offset, 1/n 
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Fig. 2.— Normalized frequency of equilibrium, i.e. of vanishing secular tidal torque, as a 
function of tidal wave number A. The asymptotic value 6e 2 is shown with the horizontal line. 



Table 1: Parameters of the tidal model for Titan as a tentative semiliquid satellite. 


Name Description Units Values 

£ .moment of inertia coefficient 0.34 

R .radius of planet m 2.576 • 10 6 

M 2 mass of the perturbed body (satellite) kg 1.3452 • 10 23 

Mi mass of the perturbing body (planet) kg 5.68 • 10 26 

a semimajor axis m 1.2219 • 10 9 

n .mean motion, i.e. 27r/P or b d _1 0.394 

e .orbital eccentricity 0.0288 

(B — A)/C .triaxiality 1.3 • 10~ 4 

Q .gravitational constant m 3 kg -1 yr -2 66468 

tm Maxwell time yr var. 

ta Andrade time yr oo 

/i .unrelaxed rigidity modulus Pa 3 • 10 9 

a .the Andrade parameter 0.2 





















be interpreted as Titan’s effective Maxwell time being longer than that value. However, other 
circumstances of purely dynamical, rather than rheological, character can prevent the capture 
of inviscid bodies into pseudosynchronism. 


4. Second condition of pseudosynchronous equilibrium 


Capture into spin-orbit resonances, including the 1:1 (synchronous) resonance, is a prob¬ 
abilistic process, conditioned by the balance of r estorat i ve and dam ping f orces similar to a 
nonlinear, driven rotating pendulum with friction ( iGoldreich fc Peale 1966, 1968j). The subt le 
conditions of a passage through spin-orbit resonances were numerically explored by flMakarov 


2012|). In general, at each resonance 9/n = 1:2, 1:1, 3:2, 2:1, ..., the resonant term of secular 


torque will be biased (or offset from zero) by the contribution of the other secular (q — j) 
terms. For example, the bias at the 1:1 resonance in secular torque is positive (see Fig. Q] 
right), mostly defined by the left-hand shoulder of the 3:2 secular term, for moderate eccen¬ 
tricities. Physically, the positive bias is generated by the fact that the perturber moves with a 
higher than n angular velocity during one half of the orbital period, and the bodies are closer 
to each other during that part of the orbit than the other half. The emerging positive, or 
acce lerating in the p rogra de sense torque, tends to spin up the synchronized perturbed body 
flMurrav fe Dermottlll9991) . On the contrary, the bias at the 3:2 resonance is usually nega¬ 
tive (decelerating) at moderate eccentricities, mostly defined by the right-hand shoulder of the 
1:1 term. It is in fact this bias that generates the possibility of non-resonant equilibrium, or 
pseudosynchronism. The bias can well be greater in magnitude than the half-amplitude of a 
resonant kink, in which case there is no equilibrium of secular torque at all (however, capture 
into that resonance may still be possible). Obviously, the bias acts as a driving force, trying to 
drag the rotating body through a resonance, whereas the near-linear inclined segment between 
the peaks acts as a damping force, or friction, working to get rid of the excess kinetic energy 
and to capture the body into the resonance. 


Returning to the useful analogy with a driven rotating pendulum, the mechanics of a 
resonance encounter can be visualized as a gradual slowing of rotation due to the secular tidal 
torque. With each cycle, the minimum angular velocity at the top of the circular path becomes 
smaller, until the pendulum stops close to the top point, no longer able to reach it. It reverses 
the motion, suddenly switching from circulation to oscillation. The first full-amplitude swing 
in the counter direction determines the fate of the pendulum. If the work done by the friction 
force is greater than the work done by the driving torque, the pendulum will not be able to 
come over the top in the opposite direction either, and it will continue to oscillate around the 
point of equilibrium with an ever decreasing amplitude due to the damping. In the numerical 
simulation shown in Fig. [3] (left), this capture into the synchronous resonance takes place 
at about 137000 days. On the other hand, when the driving torque is strong relative to the 
friction, the pendulum will pass over the top and continue to circulate in the opposite direction. 
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Fig. 3.— Numerical integration of Titan’s spin rate approaching the 1:1 spin-orbit resonance 
for two different values of orbital eccentricity. Left: capture into synchronous rotation with 
e = 0.045. Right: capture into pseudosynchronous rotation with e = 0.046. 


Capture into a pseudosynchronous equilibrium state may occur when the secular tidal 
torque averaged over one circulation period, vanishes due to the offset in relative frequency 
9/n — 1 (Fig. [TJ right) while the minimum velocity at the top of the trajectory is still greater 
than the resonant value. In other words, the pendulum still rotates in the prograde direction 
but the secular torque is balanced out. This situation is depicted in Fig. [3] (right) for the same 
set of model parameters. The mean velocity 6 remains above 1 n indefinitely long, because the 
lowest velocity over a free libration cycle can not reach the resonant value. Should the resonant 
value be reached, an immediate transition into the synchronous state happens. 


The emerging condition of pseudosynchronous equilibrium can be q uantified us ing the pre¬ 
resonant trajectory for free librations. Let us introduce, following flGoldreich fe Pealdll968l) . 


7 = 29-2M, 


(14) 


which is the sidereal libration angle multiplied by 2. Large-amplitude free librations of bodies 
with a permanent figure are driven by the triaxial torque (Eq. ED- For the case of small 
eccentricity, which is valid for the majority of close exoplanets detected so far, as well as for 
the majority of Solar system satellites, only the dominating term from Eq. [4] can be retained 
leading to 


7 ~ — 3 


B-A 

C 


n 2 G 200 


(e) sin 7 . 


(15) 


Note that this approximation can only be used for e < 0.1. From the consideration of energy 
conservation 


7 2 = 7o + 2 e 2 n 2 (l + cosy), 


(16) 


where e = y/3G 2 oo( e )(-E> — A)/C , and y 0 is the minimum libration velocity at 7 = 7 r. Equilib¬ 
rium is achieved when the kinetic energy averaged over one libration period is constant, which 
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requires the work done by the secular tidal torque to be zero, i.e., 


r* 27 T 


^TIDE d r ]f — 0. 


( 17 ) 


Per condition of pseudosynchronism number one (Section [3]), the tidal wave number of the 
resonant mode should be quite small, in which case the kvalitet in Eq. |TT] becomes a linear 
function of v 

K c {2,v)~±A 2 \, (18) 

and the secular tidal acceleration can be approximated quite well by a linear function, 


#TIDE OC — (7 — 7pseudo) • 


(19) 


Combining Eqs. [THl 03 and |TT7| obtains 


Tpseudo 


1 

2tt 

2 

7T 



(do + 2 e 2 n 2 (l + cos 7)) 1/2 d'y 


0 + 4e 2 n 2 £ 


/ 4e 2 n 2 V 

\7o + 4e 2 n 2 / 


( 20 ) 


where £ stands for the complete elliptic integral. Because of the character of the secular 
tidal torque (Fig. [U right), the planet will decelerate if the relative root-mean-square (RMS) 
velocity over one libration period is higher than 7 pS eudo- By the same argument, the planet will 
accelerate if its RMS velocity is below the point of pseudosynchronism, unless it is captured into 
the 1:1 resonance. To avoid the resonance, the minimum velocity should be positive, 70 > 0. 
In the marginal case 70 = 0, the period of free libration becomes infinitely long, because the 
pendulum stops in the unstable equilibrium 7 = 7 = 0 and 7 = 7r. The time-average spin rate, 
therefore, nullifies. Remarkably, the time-average kinetic energy or the RMS velocity do not 
nullify, in that as the minimum velocity approaches 0, 70 —> 0, the RMS velocity approaches a 
finite positive number, namely, 

4 

bpseudo t £ Tl. (21) 

TT 

Introduce a nominal eccentricity e cr j t such that the secular component of tidal torque vanishes 
at a frequency 6 e 2 rit n. As discussed in Section |3l 0 < e cr i t < e. By this definition, the 
relative spin rate 0 pseu do — n at pseudosynchronism should be equal to 6 e 2 rit n, where e cr i t — > e 
asymptotically for A —> 0. Taking into account Eq. HU we finally arrive at the second necessary 
condition of pseudosynchronus equilibrium: 


e > 


\fB 


Vsc 


A ) 


1/2 


71 


( 22 ) 


A similar condition, only different by a constant factor \/l9/2, was derived bv lGoldre ichl ( 19661) 
for the rough constant-Q tidal model. 











11 


Calculations for a model of Titan 


Titan is one o 


the la rgest sate llites in the Solar system, and a good representation of the 


class of icy moons (jHussm ann et ah 2010]). Objects like Titan can be common in exoplanet sys¬ 


tems at relatively larger separations from their hosts, bridging the gap between ice-dominated 
and terrestrial (silicate) objects. A series of numerical integrations was performed for this study 
within the uniform body approximation for a set of parameters listed in Table 1. The degre e 
of triaxiality (B — A)/C was deduced from the numbers provided by Bills & Nimmol (120 080 . 
The Andrade time was set equal to infinity, thus nullifying the effects of inelastic creep of the 
Andrade type. For most simulations, the Maxwell time was assumed to be a rather short 0.1 
days, which, for an effective shear modulus of 3 • 10 9 Pa (suggested by F. Nimmo, priv. comm. 
2013), implies a rather low effective viscosity of the dissipating layer of 2.6 • 10 13 Pa-s. This 
may still be a plausible value for a loosely bound mixture of rock and ice, as well as partially 
molten mantles. A short Maxwell time is required to realize the pseudosynchronous equilibrium 
according to the first condition (Section E} , i.e., tm v < 0.1. The asymptotic RMS velocity of 
a pseudosynchronously rotating body is 0 pseudo = (1 + 6e 2 ) n, hence the corresponding tidal 
frequency in the main mode is u — 12 e 2 n, which amounts to 0.004 d” 1 for the present day 
eccentricity, and to 0.009 d _1 for e cr i t (Section 0]). Therefore, for the Titan model in question, 
the first condition of pseudosynchronism is fulfilled for Maxwell times shorter or approximately 
equal to 10 days. We note that the first condition can be re-written for the effective viscosity 
to be 

0.1 fi 


v ;$ 


i 2e * n < 23 > 

Interestingly, for a fixed finite viscosity, the eccentricity can not be very large lest the pseu¬ 
dosynchronism equilibrium collapses by the first condition. Therefore, there is a certain interval 
of effective viscosity for pseudosynchronism to work. This interval vanishes at a characteristic 
viscosity 

O-lV^TTA 4 n A a r l l 


^crit — 


12 


0.045 


y/{B — A)/Cn 


(24) 


A)/Cn 

so that for larger values of viscosity, this type of equilibrium is practically impossible. For 
our Titan model, the critical viscosity turns out to be 2.6 • 10 15 Pa-s. This is much smaller 
than the typical effective viscosity of cold terrestrial mantles (??Earth ~ 1-3 • 10 21 Pa-s). On the 
other hand, due to the Arrhenius-type dependence of viscosity on temperature (e.g., Doremus 
20021 ). the average viscosity drops to ~ 1.3 • 10 15 Pa-s at the bottom of the Earth’s mantle 
(temperature 4000 K, pressure 136 GP af) close to the liquidus, ac cording to both the models 
presented by Shoii fe Kurita ( 2014 ) and Stamenkovic et al. ( 2012li. Similarly small pre-solidus 
viscosity is reached in the models considered by iHenn ing & Hurfordl (2 0141 ). With a fixed /i, 
planets and satellites with smaller triaxial deformation or more distant from their hosts have 
more room for a locally stable non-synchronous state, in accordance w ith Eq. [24j The triaxial 
deformation of smaller Saturnian satellites is so large OTi scareno et al.ll2009f ). that it probably 
eliminates the possibility of pseudosynchronism. 


Several numerical simulations of spin evolution for the Titan model were performed to 
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verify the second condition of pseudosynchronism (Section H]). The input parameters were 
kept as in Table [0 except for the eccentricity, which varied between 0.02 and 0.06. Each 
integration solves numerically an ordinary differential equation of second order (Eq. (2 J) for the 
angular acceleration including both the triaxial (Eq. [3]) and tidal (Eq. 0 torques. The initial 
conditions were set to 0(0) = 0, 0(0) = 1.4 n, so that the moon would start far enough from 
and above the 1:1 resonance to allow for a free, dissipative descent into an equilibrium. The 
standard Runge-Kutta integration algorithm in Mathematica was employed with an integration 
step of 0.1 days and an accuracy goal of 9 significant digits. The tidal torque included both 
the secular and oscillating terms (although the latter are probably negligible in the presence 
of much larger free librations), and the summation of tidal modes was performed over indices 
j = —1, 0,..., 4, as per Eq. [51 which is quite sufficient for such small eccentricities. The length 
of integration was set at 200000 days, which would assure a capture into one or the other 
equilibrium. 

Fig. [3] shows the resulting relative spin rate for two fixed values of eccentricity, e = 
0.045 (left) and e = 0.046 (right). This small change in the orbital eccentricity results in a 
dramatically different outcome of the encounter with the dominating resonance. In the former 
case, the spin rate at the minimum of a free libration cycle reaches the resonance 1 n, followed 
by immediate capture into that resonance. In the latter case, the point of equilibrium torque 
rotation is far enough from 1 n, the minimum spin rate can never quite reach the resonance, 
and the body remains in a super-synchronous rotation. According to Eq. [22J, the theoretical 
threshold eccentricity enabling pseudosynchronism is 0.0449. These computations show that 
the actual threshold eccentricity is between 0.045 and 0.046. The cause for this slight difference 
is in the asymptotic character of the assumed equilibrium velocity (1 + 6e 2 )n (Fig. [2]). 


3. XKT 12 ; 
2.x1(T 12 

7 1.x 1(T 12 

L. 

0 


5 

Id 


-1.x 10" 12 ;* 
-2.x10"' 
-3.x1(T 12 k 


rate of rotation 61 n 


Fig. 4.— Rate of secular eccentricity change for the Titan model (Table 1) with two fixed 
values of eccentricity: e = 0.1 (outer curve) and e = 0.0288, or the current value (inner curve). 
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6. Escapes from equilibrium 


Our considerations have been thus far decoupled from the orbital evolution processes 
caused by the secular components of the polar tidal torque. A comprehensive analysis of the 
problem should, in principle, include the total energy and angular momentum of the system 
compounded of the rotational and orbital motion of the bodies. The classic Kaula’s approach, 
on the other hand, is based on the non-inertial reference frame related to the equator of 
the perturbed body, and is not amenable to a reformulation in a barycentric non-rotating 
reference frame. Some justification to this approach is provided by the long-known fact that 
the characteristic time of s pin rate ev o lution is typically much shorte r than the time of orbit 
evolution driven by tides flZahnl 119771; lAuclair-Desrotour et al.l 120141 ). On the time scales of 
the order of 10 5 years (Fig. [3D, it is completely legitimate to ignore the secular changes in the 
orbital parameters. On the time-scales of a stellar age, however, these changes should be taken 
into account. 


The orbital eccentricity, in particular, is a crucial parameter involved in the conditions 
of pseudo synchronism . Using Kaula’s frequency decomp osition of the perturbing potential 
approach (jKaulalll964l; lLambeckljl975l; iRemus et al.ll2012n . we have for the rate of eccentricity 
from the dominating l = 2 ternj 2 


x 

x 


/de\ nVl-e 2 2 (2—m)!/ 0 X 

W ~ (Mi+M 2 )a 5 e 2^m=0 (2+m)! ° 0m 


Ep=o EjEco G>,(e) [\/l - e 2 (2 - 2p + q) - (2 - 2p)] 


M 2 R\ i&*(ii)*i(2, 4ll„) Sign(dlk) + MiRl F} mp (i 2 )K 2 (2, Sign( w «„) 


(25) 


where K { and I\ 2 are the tidal quality functions from Eq. [S] for the primary and the secondary 
components (Saturn and Titan, in our case), respectively, i\ and i 2 are the obliquities of the 
orbital plane on the respective equators, and R\ and R 2 are the radii. Ignoring the term 
corresponding to the dissipation in Saturn, which is probably smaljf], this equation implies a 
negative rate of eccentricity everywhere in the neighborhood of the 1:1 resonance for a range of 
past eccentricities, see Fig. [4j Thus, a Titan in pseudosynchronous or synchronous equilibrium 
should be slowly damping the eccentricity, barring third-body interactions. The upper bound 
equilibrium frequency, being proportional to e 2 , will decline faster, squeezing the range of 
pseudosynchronism. Eventually, the critical eccentricity is achieved and the body falls into the 
1:1 resonance. 


This process was confirmed by numerical simulations. Fig. [5] depicts the results of a 
numerical simulation similar to the ones in Fig. [3] but for a longer time, and with a variable 
eccentricity. The initial parameters were set to 0(0) = 0, 0(0) = 1.03 n, while a slowly declining 


2 For comparison, see similar equation derived for the CTL model in 


dMignardlll98c[ll98lh . 


See, however, a small Q for Saturn in (ILainev et al . [2012j) estimated from astrometric observations of its 


satellites. 
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eccentricity was modeled as 0.06 x cos(3 • 10~ 6 1), bringing it down to approximately 0.04 by the 
end of the integration. If the initial eccentricity were kept fixed during the integration span, 
the model body would descent into the pseudosynchronous state and remain there indefinitely 
long staying clear of the resonance. But the declining eccentricity pushed the RMS velocity 
down, the minimum free libration velocity reached In, and immediate capture into the 1:1 
resonance followed. 


7. Rotationally symmetric planets 

Thus far, we have considered celestial bodies which have nonzero elongation in the distri¬ 
bution of mass and in shape, when in the sequence of principal moments of inertia A < B < C, 
the smallest two are unequal, B > A. This class of objects is of special practical interest, as it 
presumably includes most Earth-sized planets and smaller satellites. Only the greater planets 
of the Solar system do not fall into this category. To a first-order approximation, moons with 
deep oceans covered by a rigid surface crust can also be described by our model, as long as 
the crust is not perfectly symmetric and is intact as a wholeB Other plausible configurations 
include planets with massive surface oceans but rigid cores, the latter likely to have nonzero 
triaxial deformation. Finally, close-in terrestrial planets can be molten to the surface by tidal 
dissipation or stellar irradiation, in which case they can not maintain a permanent figure. A 
bulge can only be raised by external perturbations, and should be dynamic and transient in 
nature. It is not obvious if the viscoelastic model in question is at all valid for such bodies. 
However, the commonly employed CTL model is mathematically similar to the Maxwell model 
in the linear regime. Therefore, the case of A = B may at least be of theoretical interest, 


4 A more accurate and consistent treatment of bodies with subsurface oceans requires more complex modeling 
including tidal wave propagation and boundary friction, best handled by numerical methods. 



0 50000 100000 150000 200000 250000 300000 

time days 


Fig. 5.— Numerically simulated transition of Titan’s spin rate from pseudosynchronous equi¬ 
librium to the 1:1 (synchronous) resonance caused by a gradual decline of orbital eccentricity 
from 0.06 at t — 0 to 0.042 at t = 300000 days. 
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delineating the boundaries of our knowledge within the existing analytical models. 

Obviously, the second condition of pseudosynchronism (Eq. [22} is automatically satisfied. 
At first glance, the first condition also appears defunct, because any finite separation of the 
point of equilibrium secular torque from exact resonance should guarantee capture into this 
equilibrium in the absence of oscillations. Indeed, there is no triaxial torque, thus, there is 
no free libration. A more careful analysis, supported by numerical experiments, shows that 
other, much smaller librations are still present. These are always harmonics of the orbital 
frequency caused by the periodic terms of the tidal torque, Eq. 0 We observe that the mixed 
terms q ^ j give both sine- and cosine-harmonics of perturbing polar torque, and hence, similar 
variations of the spin rate and libration angle. For a given pair q, j, the difference between these 
harmonics is in the kvalitet functions K c and K s . Recalling that GW- 1 ) = —1/2 e + 0(e 3 ), 
G 200 = 1 — 5/2 e 2 + 0(e 4 ), G 201 = 7/2 e + 0(e 3 ), and all other Kaula functions are 0(e 2 ), or 
higher orders in e, our analysis can be safely limited to the combinations of 0 with -1 or 1 in 
the (q,j) indices, for a small eccentricity. Furthermore, as long as the spin rate is close to the 
1:1 resonance, the largest by magnitude terms among the cos-harmonics will be the ones with 
q — ±1 and j = 0, because the q — 0 term has K c (2, 0) = 0. The kvalitet for these two terms 
is the same, equal to K c (2,n) (Eq. [IT]). For the case M 2 <C M u this leads to an approximate 
equation for the combined cos-harmonic of tidal forced libration 


0 ~ 
w cos 


a M 1 2 

6 ——n 

iM 2 


R 


K c ( 2, n) e cos Ad. 


(26) 


The corresponding forced libration in angular velocity is 9 — n oc sin Ad, and in angle 6 — Ad oc 
— cosAd. Apart from being numerically very small, this perturbation can not violate the 
pseudosynchronous equilibrium, because it is out of phase with the orbital acceleration of 
the perturber. Indeed, in the absence of a triaxial force, any restorative force should still be 
proportional to + sin Ad to lock the body into synchronous rotation. 


The ' sin-harmonics of oscillating tidal torque emerge in the general Darwin-Kaula expan¬ 
sion (lEfroimskvll2012al) . but, to my knowledge, have not been exploited in the literature. They 
are proportional to the kvalitet functions K s (2,is) (Section [ 2 ] Eq. [5]), which are, explicitly, 


_ 3 A 2 ut m ^R(u) + 3ft(z/) 2 + A(R) 2 
2 (!R(z/) + A 2 vtm) 2 + A 2 (^) 


(27) 


For negligible Andrade creep contribution (ta = 00 ), this further simplifies to 


K s ( 2,u) 


3 Al 2 A 2 + A 2 + 1 
2 A 2 (l + A 2 ) 2 + 1' 


(28) 


where A = utm is the tidal wave number. The maximum value of K s is reached at u — 0, 
K s ( 2,0) = 3/2, whereas K c (2, u 2200 ) ~ 3/2 .A 2 ^2200 for A^^ooTm <C 1. This indicates that 
the sin-harmonics may in fact dominate the spectrum of tidal libration. There is, however, 
another circumstance that may prevent this. For very low effective viscosities and sufficiently 
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wide orbits, when utm <C 1, the kvalitet is a slowly declining function of frequency, 

K s {2. u) — — (l — > 4 . 2(1 + >4. 2 )A 2 + 0(A 4 )) , 


( 29 ) 


so that K s (2,n) ~ K s (2,0). In this case, the pairs of symmetric terms q,j = 0,1 and 1,0, and 
0 , —1 and —1,0, almost cancel each other, being opposite in sign, and differing by the factor 
K s ( 2, v) in magnitude. This happens to our model Titan with a nominal tm = 0.1 d considered 
in the previous simulations. For the sake of completeness, we should consider another specific 
scenario when A 2 nrM 3> 1. Then, K s (2,ri) ~ 3/2(1 + >4 2 ) and, limiting the expansion in e to 
the first power, 


6 • 

w sin ^ 


R Ml 2 
b—~r n 
£M 2 


R 


A 2 


e sin A4. 


(30) 


1 + 4l 2 

Due to the positive sign of this periodic force, it acts in a similar fashion to the triaxial torque 
as a restorative force. A stable pseudosynchronous equilibrium is not possible if the offset of 
the equilibrium torque in frequency is smaller that the amplitude of the forced libration caused 
by this periodic tidal torque, which leads to 


Mi (R\‘ A 2 

(M2\aJ 1+-V 


(31) 


and the only condition of pseudosynchronism is e > e cr ; t . For Titan, the calculated value of e cr ; t 
is a rather low 9 • 10~ 5 . Such small numbers are expected for practically all known star-planet 
and planet-moon systems because R a. In conclusion, we identified one area of parameters 
for rotationally symmetric bodies (>4 2 nTM 1) where pseudosynchronous equilibrium can be 

precluded by the oscillating components of the tidal torque, but the corresponding threshold 
eccentricity is so low that such cases should hardly exist. 


8. Summary 


The rate of energy dissipated in a uniform deformable body by tides in the dominating 
degree / = 2 is 


f dE 
v dt ' 


QMI& 2 (2-m)! 

/ , 20 , ™m( 2 - " 0m ) 


a u 


(2 + m)\ 


m =0 
2 +00 

-EE [-^2mp(^) Cr2pq{^)\ ^ / 2mpq-^-c{^ / 2‘ 

p =0 q=—oo 


rnpq ) 


(32) 


where i is the equator obliquity relative to the orbit plane, Fv m. P (i) a re the inclination functions, 
and other parameters are defined elsewhere flEfroimskv fc Makarovl l2014l). Due to the presence 
of the positively definite tidal frequency V 2 mpq , the principal resonant mode of the tidal torque 
is nullified and does not participate in the dissipation of kinetic energy. For example, the main 
component of the torque for small and moderate eccentricities is the mpq = 200 term, but it 
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vanishes in Eq. [32] if the planet is captured into the 1:1 spin-orbit resonance. The other tidal 
modes, separated from the dominant mode by an integer times n, will continue to contribute 
to dissipation. A semiliquid body caught in pseudosynchronous rotation, on the other hand, 
continues to rotate with respect to the perturber’s vector in the prograde sense, so that the 
dominant mode of tidal torque continues to generate some dissipation (limited due to the 
smallness of the tidal frequency ^ 2200 ); however, the negative torque from this mode will be 
exactly balanced, on average, by the component torques from the other modes. As a result, the 
rate of tidal dissipation is smaller at pseudosynchronism than at synchronism. The fact that 
the heating rate is minimal at the point of tidal equilibrium has been known and discussed in 
the literature ( Wisdom 2008; Levrard 2008). 


I 11 this paper, only the “principal” pseudosynchronous state has been considered explicitly, 
i.e. the one associated with the 1:1 spin-orbit resonance. This limits the analysis to the case 
of small and moderate eccentricities. The same mathematics can be easily generalized for 
larger eccentricities (e > 0 .2), where e quilibrium torque emerges at other tidal modes. Using 
the terminology by (IStorch fe Lai 201411 . only the true, or resonant, equilibria will be stable, 
whereas unstable equilibria of null secular torque can be found between the resonances. These 
have no practical significance, and they are not considered in this paper. This may be true 
only for purely Maxwell-type viscoelastic models without any secondary inelastic mechanisms 
of dissipation, however. The action of an Andrade-type dislocation and unpinning mecha nism, 
for example, coupled with a characteristic low-frequency cutoff (IKarato fe Spetzler 19 901). may 
concievably lead to locally stable equilibria between the spin-orbit resonances. Observations 
of the Lunar tidal phase lag with GRAIL indicate the existence of a secondary peak in th e 
general kvalitet factor between the frequencies 0.075 n and In ( Williams et al.l 2015a ). 


This secondary peak can be heuristically fitted with a Debye peak model ([Williams et ah 


2 015b c). but the Andrade creep model with a ta — 0.45 yr can also reproduce this feature. 
If the Moon’s eccentricity were higher, the secondary peak of the main resonance could be 
shifted toward positive values of secular torque by the left shoulder of the 3:2 resonance in 
such a way that a zero-torque crossing with a locally negative slope takes place between these 
two resonances. This should make non-resonant rotational equilibria possible in principle. 
A more detailed analysis supported by numerical simulations is required to investigate this 
curious possibility. 


We have demonstrated that pseudosynchronous rotation is only possible under specific, 
restrictive conditions. As long as the assumption of a homogeneous body is applicable, it 
can only happen to rather inviscid, semiliquid or semi-molten bodies with a characteristic 
tidal wave number A = tm n -Cl. If this condition is not satisfied, the offset of equilibrium 
secular torque collapses and complete synchronization is the only stable state. Secondly, the 
orbital eccentricity should be large enough to accommodate the range of free libration to avoid 
overshooting the resonance at the minimum velocity, per Eq. [22 ] The threshold value depends 
on the degree of elongation (B — A)/C, hence, more rotationally symmetric bodies have better 
chances to be pseudosynchronous than their more deformed counterparts. On the other hand, 
eccentricity may be limited by the first condition for fixed viscosity and shear modulus, see. 
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Eq. [24j The combination of these opposing criteria implies that pseudosynchronous equilibria 
are practically impossible for large satellites like Titan if the effective viscosity is greater than 
« 10 15 Pa-s. 

One possibly observable aspect of a stable pseudosynchronous state is that the large- 
amplitude, long-period free librations of angular velocity around the equilibrium value 7 pse udo 
are eternal, i.e., not diminishing in amplitude with time (cf. Fig. [3l right). Physically, this 
fact is readily explained by the null work produced by the secular tidal torque in that state. 
For the free libration to be damped, the work of secular torque needs to be negative. This 
happens when the periodic torque and velocity have opposite phase, whereas in our case, the 
torque changes its sign over one period of libration but the velocity 7 remains positive. On the 
other hand, the tidal bulge constantly moves across the surface, the body is being deformed 
and stressed, and the energy dissipation, being minimized at pseudosynchronism, still goes on. 
Obviously, the energy spent on tidal heating is not taken from the kinetic energy of the rotating 
body, but directly from the orbital energy. Inevitably, the orbit should gradually shrink, unless 
it is supported by a fast rotating perturber. 

Being only marginally stable, a pseudosynchronous body can spontaneously leave this 
state and fall into the 1:1 resonance if the eccentricity falls below the critical value. The 
tidally driven rate of eccentricity is negative in the vicinity of the 1:1 resonance for small 
e. Most tidally interacting planets and satellites should therefore be eventually captured in 
that resonance, unless the eccentricity is boosted by the tides in the primary or other bodies. 
This may be the reason we do not find any pseudosychronously rotating satellites in the Solar 
system. Titan, with its current value of eccentricity, could not be pseudosynchronous, but if 
its orbit was more elongated in the past, it could be pseudosynchronous for a limited duration. 
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